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ABSTRACT 

Asteroseismology of bright stars with well-determined properties from parallax measurements and interferom¬ 
etry can yield precise stellar ages and meaningful constraints on the composition. We substantiate this claim 
with an updated asteroseismic analysis of the solar-analog binary system 16 Cyg A & B using the complete 30- 
month data sets from the Kepler space telescope. An analysis with the Asteroseismic Modeling Portal (AMP), 
using all of the available constraints to model each star independently, yields the same age (t = 7.0 ±0.3 Gyr) 
and composition (Z = 0.021 ±0.002, Y, = 0.25 ±0.01) for both stars, as expected for a binary system. We 
quantify the accuracy of the derived stellar properties by conducting a similar analysis of a Kepler-like data 
set for the Sun, and we investigate how the reliability of asteroseismic inference changes when fewer observa¬ 
tional constraints are available or when different fitting methods are employed. We find that our estimates of 
the initial helium mass fraction are probably biased low by 0.02-0.03 from neglecting diffusion and settling 
of heavy elements, and we identify changes to our fitting method as the likely source of small shifts from our 
initial results in 2012. We conclude that in the best cases reliable stellar properties can be determined from 
asteroseismic analysis even without independent constraints on the radius and luminosity. 

Keywords: stars: individual (HD 186408, HD 186427)—stars: interiors—stars: oscillations—stars: solar-type 


1. INTRODUCTION 

Asteroseismology is emerging as a powerful technique to 
determine the ages and other properties of field stars (Met¬ 
calfe et al. 2014), including the parent stars of planetary sys¬ 
tems (Silva Aguirre et al. 2015). During the nominal Kepler 
mission, the data analysis and modeling methods for aster¬ 
oseismic inference matured substantially, and the techniques 
are now being applied to brighter stars and members of well- 
characterized clusters for the K2 mission (Howell et al. 2014; 
Chaplin et al. 2015). These studies will ensure that automated 
methods to determine stellar properties are as reliable as pos¬ 
sible prior to the collection of asteroseismic data for bright 
stars all around the sky with the Transiting Exoplanet Survey 
Satellite (TESS, Ricker et al. 2014). 

The old solar analogs 16 Cyg A & B are among the bright¬ 
est stars (V ~ 6) within the original Kepler field of view. Pre¬ 
cise photometric measurements were obtained every 58.85 
seconds almost continuously for 2.5 years, between 2010 
September 23 and 2013 April 8 (quarters Q7-Q16). Met¬ 
calfe et al. (2012) conducted an initial asteroseismic analy¬ 
sis of 16 Cyg A & B using the first 3 months of Kepler data, 
identifying more than 40 oscillation modes in each star. The 
optimal models, derived by fitting the observations of each 
star independently, had the same age and initial composition 
within the uncertainties—as expected for the components of a 
binary system. These initial results bolstered our confidence 
in the reliability of asteroseismic techniques. 

Since the original study by Metcalfe et al. (2012) the time 
series for 16 Cyg A & B has been extended by an order 
of magnitude, and radius constraints from interferometry are 
now available (White et al. 2013). The complete Kepler data 
sets have already been used to constrain the bulk helium abun¬ 
dances (Verma et al. 2014), and to determine rotation rates 
and inclinations (Davies et al. 2015) for both components. 


Using these new observational constraints, and applying the 
improved modeling methods developed over the past several 
years, the aim of this Letter is to quantify the precision and 
accuracy of asteroseismic inferences for the bright targets that 
will be observed by future missions such as TESS and PLATO 
(Rauer et al. 2014). 

In Section 2 we outline the previously published observa¬ 
tions that we adopt for our analysis, including a photometric 
data set for the Sun degraded to Kepler-like precision. We 
summarize recent updates to our methods in Section 3, and 
we describe a new set of experiments using the Asteroseismic 
Modeling Portal (AMP, Woitaszek et al. 2009; Metcalfe et al. 
2009) to probe the influence of different observational con¬ 
straints and fitting methods. Finally, in Section 4 we discuss 
the precision and accuracy of the derived stellar properties for 
16 Cyg A & B, using the solar results to identify any biases. 

2. OBSERVATIONAL CONSTRAINTS 

To update our asteroseismic analysis, we adopted the sets 
of frequencies published by Davies et al. (2015), including 
54 and 56 oscillation modes for 16 Cyg A & B respectively. 
These frequency sets were extracted from the power spectra 
of 928 days of short-cadence observations with a duty cycle of 
90.5 percent. Each set includes 15 consecutive orders for the 
radial (/ = 0), dipole {l = 1) and quadrupole (/ = 2) modes, as 
well as 9 and 11 octupole (/ = 3) modes for 16 Cyg A & B 
respectively. Compared to the 3 months of data presented 
in Metcalfe et al. (2012), the longer time series analyzed by 
Davies et al. (2015) typically allowed the detection of 1-2 ad¬ 
ditional orders at lower and higher frequencies and improved 
the precision by a factor of 2-4 for previously detected modes, 
consistent with the expectations of Libbrecht (1992). 

In addition to the frequencies, we adopted other obser¬ 
vational constraints from spectroscopy, Hipparcos parallaxes 
(van Leeuwen 2007), and interferometry. The spectroscopic 


2 


METCALFE ET AL. 


constraints and luminosities were similar to those used by 
Metcalfe et al. (2012): T et£A = 5825 ±50 K, [M/H] A = 0.10± 
0.09, T e f fB = 5750 ±50 K, [M/H]* = 0.05 ±0.06 (Ramirez 
et al. 2009), L A = 1.56 ± 0.05 L e and L B = 1.27 ± 0.04 L e 
(Metcalfe et al. 2012). Note that we adopted 3cr uncertain¬ 
ties on the spectroscopic metallicities to allow for potential 
systematic errors. Unlike the 2012 study, radius constraints 
are now available from the interferometric observations of 
White et al. (2013), who used the CHARA array (ten Brum- 
melaar et al. 2005) with the PAVO beam combiner (Ireland et 
al. 2008). They found linear radii R A = 1.22 ±0.02 Rq and 
Rb = 1.12 ±0.02 Rq, which were combined with an aster- 
oseismic scaling relation (Ulrich 1986) to obtain mass esti¬ 
mates M a = 1.07 ± 0.05 Mq and M B = 1.05 ± 0.04 M 0 . The 
radius estimates are independent of asteroseismology, so we 
adopt them as constraints for the modeling presented in Sec¬ 
tion 3. We note the mass estimates here only for comparison 
with our final results (see discussion in Section 4). 

To assess any biases in our modeling, we adopted a pho¬ 
tometric data set for the Sun that was used by Davies et al. 
(2015) for validation of their analysis methods. The solar 
data were obtained from the red channel of the VIRGO instru¬ 
ment (Frohlich et al. 1995) on the SoHO satellite (Domingo 
et al. 1995) with white noise added to approximate the Kepler 
data for 16 Cyg A & B. The solar power spectrum was ana¬ 
lyzed using the same methods that were used for 16 Cyg, and 
the extracted frequencies were chosen to include the same set 
of modes (relative to the frequency of maximum power) de¬ 
tected in 16 Cyg A. The resulting solar frequency errors are 
comparable to those obtained for 16 Cyg A, and all other ob¬ 
servational constraints for the Sun were fixed at the known 
values with the same fractional errors as 16 Cyg A: T e g 0 = 
5777 ±50 K, [M/H]q = 0.00 ±0.09, L 0 = 1.000 ±0.032, and 
Rq = 1.000 ± 0.016. 

3. ASTEROSEISMIC MODELING 

We used the observational constraints described in Sec¬ 
tion 2 to obtain the optimal stellar properties with AMP from 
various data sets and fitting methods. The AMP science code 
uses a parallel genetic algorithm (Metcalfe & Charbonneau 
2003) to optimize the match between models produced by 
the Aarhus stellar evolution and pulsation codes (Christensen- 
Dalsgaard 2008a,b) and a given set of observations. The mod¬ 
els are configured to use the OPAL 2005 equation of state 
(Rogers & Nayfonov 2002) with opacities from OPAL (Igle- 
sias & Rogers 1996) and Ferguson et al. (2005), and nu¬ 
clear reaction rates from the NACRE collaboration (Angulo 
et al. 1999, 2005). Convection is treated using mixing-length 
theory (Bohm-Vitense 1958) without overshoot, while diffu¬ 
sion and settling of helium is treated with the prescription of 
Michaud & Proffitt (1993). The AMP software has been in 
development since 2004, and below we briefly outline its his¬ 
tory to place the current modeling approach in context. 

3.1. Updated physics and methods 

The AMP 1.0 software (Metcalfe et al. 2009) reproduced 
the solar properties from Sun-as-a-star data, but we had more 
difficulty fitting the early asteroseismic observations from Ke¬ 
pler (Metcalfe et al. 2010). For AMP 1.1 (Metcalfe et al. 
2012; Mathur et al. 2012) we used the statistical errors on 
each oscillation frequency to assign weights, rather than us¬ 
ing a combination of statistical and systematic errors that 
gave lower weight to higher frequencies where so-called “sur¬ 
face effects” dominated. We also split our % 2 quality metric 


into separate components for seismic and non-seismic observ¬ 
ables, to prevent the many oscillation frequencies from over¬ 
whelming the relatively few and less precise spectroscopic 
constraints. For AMP 1.2 (Metcalfe et al. 2014) we updated 
the model physics to use the NACRE reaction rates (Angulo 
et al. 1999) instead of those from Bahcall et al. (1995) and the 
low-temperature opacities of Ferguson et al. (2005) instead 
of Alexander & Ferguson (1994). We also began fitting two 
sets of frequency ratios row and /yo (Roxburgh & Vorontsov 
2003) in addition to the frequencies themselves, combining 
these three sets of observables with a set of spectroscopic and 
other constraints into an average \ 2 quality metric. 

We have made several important updates to AMP for the 
results presented here. First, we implemented the revised rate 
for the 14 V +p reaction determined by the NACRE collabora¬ 
tion (Angulo et al. 2005). Second, we now use only the fre¬ 
quency ratios row and ro 2 (along with r 13 for 16 Cyg A & B 
where the / = 3 modes are detected) rather than the individual 
frequencies for modeling. This allows us to avoid potential bi¬ 
ases from applying an empirical correction for surface effects 
(e.g. Kjeldsen et al. 2008). The frequency ratios are insensi¬ 
tive to the near-surface layers by design, so they allow a direct 
comparison between observations and models without resort¬ 
ing to any ad hoc frequency corrections. The frequency ra¬ 
tios are also insensitive to line-of-sight Doppler velocity shifts 
(see Davies et al. 2014). Finally, we combine all of the obser¬ 
vational constraints into a single y 2 quality metric to facili¬ 
tate the statistical interpretation of the results. The average 
X 2 metric used in AMP 1.2 was designed to yield the optimal 
trade-off between complementary sets of constraints, but it 
complicated the determination of reliable uncertainties on the 
inferred stellar properties. We calculated the uncertainties in 
the same manner as described in Metcalfe et al. (2014), but us¬ 
ing the single y 2 value instead of an average y 2 to determine 
the likelihood of each sampled model. The AMP 1.2 fitting 
method can optionally still be used, and we have implemented 
the scaled solar surface correction proposed by Christensen- 
Dalsgaard (2012) to help minimize any resulting biases. We 
demonstrate below that our results are largely insensitive to 
the particular choice of fitting method. 

3.2. Results for 16 Cyg A & B and the Sun 

We used AMP to investigate how the reliability of astero¬ 
seismic inference changes when various sets of observational 
constraints are adopted or when slightly different fitting meth¬ 
ods are employed. The optimal asteroseismic properties of 
16 Cyg A & B and the Sun for each case are listed in Table 1, 
and the variations are described below. 

The most reliable properties, obtained by applying the up¬ 
dated fitting method described above to the complete set of 
observational constraints 1 , are shown in the first case labeled 
“all constraints” in Table 1 (AMP simulations 767, 768 and 
766 for 16 Cyg A, 16 Cyg B and the Sun respectively). All 
three models match the spectroscopic and other constraints 
with a reduced y 2 < 1. The quality of the match to the astero¬ 
seismic constraints for 16 Cyg A & B is illustrated in Figure 1 
where the bottom panels show the match to the frequency ra¬ 
tios that were actually used as constraints for the modeling, 
while the top panels show the match to the frequencies after 
the scaled solar surface correction of Christensen-Dalsgaard 

1 The constraints used for each case are provided as an obs.dat file on the 
AMP website at http://amp.phys.au.dk. 
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Table 1 

Properties of the optimal models using various constraints and fitting methods 


Case 

R/R @ 

M/M© 

L/L© 

t! Gyr 

Z 

Y, 

l°g£ 

AMP a 

16 Cyg A (Kepler) 

1: all constraints. 

1.229 ±0.008 

1.08 ±0.02 

1.55 ±0.07 

7.07 ±0.26 

0.021 ±0.002 

0.25 ±0.01 

4.292 ±0.003 

767 

2: without R. 

. 1.225 ±0.008 

1.07 ±0.02 

1.53 ±0.07 

7.15 ±0.27 

0.020 ± 0.002 

0.25 ±0.01 

4.291 ±0.002 

773 

3: without R,L. 

1.225 ±0.007 

1.07 ±0.02 

1.52 ±0.07 

7.12 ±0.23 

0.021 ±0.001 

0.25 ±0.01 

4.291 ±0.002 

776 

4: fitting method. 

1.229 ±0.007 

1.08 ±0.02 

1.55 ±0.06 

7.07 ±0.37 

0.021 ±0.002 

0.25 ±0.01 

4.292 ±0.002 

770 

5: surface term. 

1.229 ±0.008 

1.08 ±0.02 

1.55 ±0.07 

7.07 ±0.34 

0.021 ±0.002 

0.25 ±0.01 

4.292 ± 0.002 

797 

6: 3 months data. 

1.229 ±0.008 

1.08 ±0.02 

1.55 ±0.07 

6.99 ±0.37 

0.024 ±0.002 

0.26 ±0.01 

4.292 ±0.003 

778 

16 Cyg B (Kepler) 

1: all constraints. 

1.116 ±0.006 

1.04 ±0.02 

1.25 ±0.05 

6.74 ±0.24 

0.022 ± 0.003 

0.26 ±0.01 

4.359 ±0.002 

768 

2: without R. 

. 1.116 ±0.006 

1.04 ±0.02 

1.25 ±0.06 

6.74 ±0.23 

0.022 ±0.002 

0.26 ±0.01 

4.359 ±0.002 

774 

3: without R,L. 

1.116 ±0.006 

1.04 ±0.02 

1.24 ±0.06 

6.79 ±0.19 

0.022 ± 0.002 

0.26 ±0.01 

4.359 ±0.002 

777 

4: fitting method. 

1.116 ±0.005 

1.04 ±0.01 

1.25 ±0.05 

6.89 ±0.28 

0.020 ± 0.002 

0.25 ±0.01 

4.359 ±0.002 

771 

5: surface term. 

. 1.116 ±0.006 

1.04 ±0.02 

1.24 ±0.05 

6.96 ±0.29 

0.020 ±0.002 

0.25 ±0.01 

4.359 ±0.002 

798 

6: 3 months data. 

1.113 ±0.006 

1.03 ±0.02 

1.25 ±0.05 

6.66 ±0.27 

0.022 ±0.001 

0.27 ±0.01 

4.358 ±0.002 

779 

Sun (SoHO/VIRGO) 

1: all constraints. 

1.003 ±0.006 

1.01 ±0.02 

0.98 ±0.04 

4.62 ±0.15 

0.018 ±0.002 

0.26 ±0.02 

4.439 ± 0.002 

766 

2: without R. 

. 1.007 ±0.006 

1.02 ±0.02 

1.02 ±0.04 

4.53 ±0.15 

0.018 ±0.002 

0.26 ±0.01 

4.440 ±0.002 

772 

3: without R,L. 

. 1.010 ±0.005 

1.03 ±0.02 

1.00 ±0.04 

4.58 ±0.14 

0.018 ±0.001 

0.25 ±0.01 

4.442 ± 0.002 

775 

4: fitting method. 

1.010 ±0.007 

1.03 ±0.02 

1.00 ±0.04 

4.52 ±0.18 

0.020 ±0.003 

0.25 ±0.01 

4.442 ± 0.002 

769 

5: surface term. 

1.003 ±0.007 

1.01 ±0.02 

1.00 ±0.03 

4.61 ±0.18 

0.018 ±0.003 

0.26 ±0.02 

4.439 ±0.002 

799 

6: 3 months data. 

. 0.996 ±0.010 

0.99 ±0.03 

1.00 ±0.04 

4.77 ±0.18 

0.017 ±0.001 

0.27 ± 0.02 

4.437 ± 0.004 

800 


a Comprehensive model output is available at http://amp.phys.au.dk/browse/simulation/### 


(2012) has been applied to the raw model frequencies. All of 
the models that result from variations to the set of observa¬ 
tional constraints and to the fitting methods yield comparable 
matches to both the asteroseismic and other constraints. The 
primary reason for exploring such variations is to quantify the 
absolute accuracy of the inferred stellar properties, and to as¬ 
sess the sensitivity of the results to details of the modeling 
strategy and data quality. 

We can test the absolute accuracy of our results by elim¬ 
inating the independent constraints on radius and luminosity 
to see whether the asteroseismic and spectroscopic constraints 
alone are sufficient to determine these properties. Interfer¬ 
ometric radius constraints are available only for relatively 
nearby stars, so the second case (labeled “without R” in Ta¬ 
ble 1, AMP simulations 773, 774 and 772) removes this con¬ 
straint from the % 2 value that guides the search. For 16 Cyg A 
and the Sun this leads to small shifts within the quoted uncer¬ 
tainties for the optimal stellar properties, while only the un¬ 
certainties change for 16 Cyg B. Luminosity constraints are 
currently available for a small fraction of stars in the Kepler 
field, at least until parallaxes are released by the Gaia mission 
(Perryman et al. 2001). Parallaxes are also needed to convert 
the angular measurements from interferometry into a linear 
radius, so our third case (labeled “without R,L” in Table 1, 
AMP simulations 776, 777 and 775) omits both the radius 
and luminosity constraints. Again there are marginal shifts in 
the optimal stellar properties, insignificant for 16 Cyg A & B 
but slightly larger than the uncertainties for the solar radius 
and mass. In general, the reliability of asteroseismic infer¬ 
ence seems to degrade more significantly for the Sun than for 
16 Cyg A & B when the radius and luminosity constraints are 
not available. 

Large samples of Kepler stars have previously been fit using 
slightly different methods to assess the quality of the match 
between the models and observations. To evaluate the impact 
of these differences on the inferred stellar properties, we fit 


the full set of constraints from the first case using two alter¬ 
native modeling strategies. For the fourth case (labeled “fit¬ 
ting method” in Table 1, AMP simulations 770, 771 and 769) 
we used the AMP 1.2 fitting method (Metcalfe et al. 2014), 
which calculates an average x 2 from four sets of observables 
including the frequencies corrected for surface effects with a 
solar-calibrated power law following Kjeldsen et al. (2008). 
Comparing the inferred stellar properties to the first case, 
16 Cyg A yields slightly different uncertainties, 16 Cyg B 
is shifted within the uncertainties, and only the Sun shows 
marginally significant shifts in radius and mass. To probe the 
possible source of these biases, for the fifth case (labeled “sur¬ 
face term” in Table 1, AMP simulations 797,798 and 799) we 
repeated the fourth case using the scaled solar surface correc¬ 
tion of Christensen-Dalsgaard (2012) to replace the Kjeldsen 
et al. (2008) prescription. The differences were insignificant 
for 16 Cyg A & B, but for the Sun the alternative surface cor¬ 
rection effectively eliminated the biases. 

Finally, we examine the influence of using the longer aster¬ 
oseismic time series by applying our updated fitting method to 
the original data sets for 16 Cyg A & B published in Metcalfe 
et al. (2012), along with a comparable 3-month solar data 
set from VIRGO. This final case (labeled “3 months data” 
in Table 1, AMP simulations 778, 779 and 800) should be 
compared to the second case “without R”, since the original 
data sets did not include the radius constraints of White et al. 
(2013). For 16 Cyg A & B these two cases yield similar stel¬ 
lar properties, indicating that differences relative to the Met¬ 
calfe et al. (2012) results can be attributed mostly to the up¬ 
dated model physics rather than the extended time series. For 
the Sun, the properties inferred from the shorter time series 
reproduce the solar radius, mass and luminosity more accu¬ 
rately. This seems to be a consequence of the fewer and lower- 
precision asteroseismic constraints, which gives the spectro¬ 
scopic and other constraints more relative weight when all of 
the observations are combined into a single % 2 quality metric. 
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16 Cyg A 


16 Cyg B 



Frequency modulo 103.4 (/r.Hz) 



Frequency modulo 116.9 (/xHz) 




Figure 1 . Echelle diagrams (top) and frequency ratios (bottom) for 16 Cyg A (left) and 16 Cyg B (right), with the observations shown as connected points 
with errors. Top panels: frequencies of the optimal models from AMP are shown using different symbols to indicate radial (circles), dipole (triangles), 
quadrupole (squares) and octupole (diamonds) oscillation modes. The scaled solar surface correction of Christensen-Dalsgaard (2012) was applied to all model 
frequencies for illustration. Bottom panels: frequency ratios roio (circles), 7*02 (squares) and r 13 (triangles) that were actually used as constraints for the modeling. 


4. DISCUSSION 

We have determined precise asteroseismic properties for the 
solar-analog binary system 16 Cyg A & B, and we have quan¬ 
tified the accuracy of the results by varying the input con¬ 
straints and fitting methods, and by applying the same tech¬ 
niques to the Sun. We obtain a radius precision of 0.6%, a 
mass precision of 2%, and an age precision of 3%. Our results 
for the Sun suggest that the systematic errors are smaller than 
or comparable to the statistical uncertainties, consistent with 
the conclusions of previous studies (e.g. Lebreton & Goupil 
2014; Silva Aguirre et al. 2015). The radii and masses of 
16 Cyg A & B are consistent with the estimates of White et 
al. (2013) from interferometry and scaling relations, but sys¬ 
tematically lower than the original AMP results in Metcalfe 
et al. (2012). These differences persist when the radius con¬ 
straint is omitted, or when the 2012 data set is analyzed with 
updated methods, suggesting that the use of frequency ratios 
as asteroseismic constraints is responsible. The luminosities 
of 16 Cyg A & B agree with the values derived by Metcalfe 
et al. (2012) from the Hipparcos parallaxes, though the best 
models for both components have luminosities slightly below 
the constraints. The solar radius, mass and luminosity are 
all recovered faithfully within the quoted precision, including 
small biases of +0.3% in radius, +1% in mass and -2% in lu¬ 
minosity. Note that results for significantly more massive or 
more evolved stars may not be as precise or accurate, even 
with comparable data quality. 


Although we fit the observations of 16 Cyg A & B indi¬ 
vidually, the results for the two stars have a common age 
(t = 7.0 ± 0.3 Gyr) and composition (Z = 0.021 ± 0.002, T; = 
0.25 ±0.01) within the uncertainties, as expected for the com¬ 
ponents of a binary system. The updated age is slightly older 
and more precise than the result presented in Metcalfe et al. 
(2012), and the updated methods in AMP now yield a consis¬ 
tent age for both components. The closest agreement between 
the two sets of models is obtained when using the AMP 1.2 
fitting method with the scaled solar surface correction of 
Christensen-Dalsgaard (2012). The comparable fits to solar 
data agree with the seismic age of the Sun (4.60 ± 0.04 Gyr) 
as determined by Houdek & Gough (2011), with the ensem¬ 
ble of models in Table 1 showing a total age spread of only 
0.25 Gyr (5%). 

As with the original results of Metcalfe et al. (2012), we 
find values for the initial helium mass fraction that are system¬ 
atically low compared to other estimates. Verma et al. (2014) 
used deviations from uniform frequency spacing in the oscil¬ 
lation modes to constrain the acoustic depths of the base of the 
convection zone (BCZ, t cz ) and the helium ionization zone. 
The amplitude of deviations due to the latter source can be 
calibrated with models to determine the current helium mass 
fraction in the stellar envelope. The weighted average of three 
different methods using two different models for calibration 
yields Y s = 0.241 ± 0.004 for 16 Cyg A and Y s = 0.244 ± 0.004 
for 16 Cyg B. In our models X, is systematically lower by 
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0.02-0.03, confirming that we are biased toward low initial 
helium. This issue had previously been identified by Gruber- 
bauer et al. (2013) from Bayesian modeling of the 2012 data 
set for 16 Cyg A & B, and our results for the Sun show the 
same problem. We speculate that the bias toward low ini¬ 
tial helium may be due to neglecting diffusion and settling of 
heavy elements, which is not stable under all conditions for 
the models used in AMP. We hope to remedy this issue with 
future development of AMP 2.0, which will use the MESA 
stellar evolution code (Paxton et al. 2011, 2013, 2015) and 
the GYRE pulsation code (Townsend & Teitler 2013). 

An additional diagnostic of our best models is whether they 
agree with the convection zone depths derived by Verma et al. 
(2014). Although unrealistic surface boundary conditions in 
the models hinder direct comparisons, Mazumdaret al. (2014) 
suggest comparing the fractional acoustic radius of the BCZ 
from observations and models. The total acoustic radius of 
the star To can be estimated from the mean frequency spac¬ 
ing of radial modes Ao using Tq ~ 1/(2Ao). The fractional 
acoustic radius of the BCZ is then I nc//To = 1 — (t cz /7o). The 
comparable quantity for a model is obtained from a numer¬ 
ical integration of the inverse sound speed from the center 
of the model to the BCZ, relative to a full integration out 
to the photosphere. A weighted average of the values of 
Tkcz/Tq from three methods using the data in Tables 1 and 
3 of Verma et al. (2014) gives 0.368 ±0.012 for 16 Cyg A 
and 0.391 ±0.028 for 16 Cyg B, while our best models yield 
0.386 and 0.381 respectively—slightly above the constraint 
for 16 Cyg A (+1.5cr), but consistent for 16 Cyg B (-0.4er). 
Future efforts should consider how to include the fractional 
acoustic radius of interior features as additional constraints 
for the modeling. 

Automated asteroseismic modeling has advanced signifi¬ 
cantly in the past few years, and we are on track to take full 
advantage of the large data sets that will emerge from future 
missions like TESS (1-12 month time-series) and PLATO (5- 
36 months per field). The analysis presented in this Letter 
demonstrates that high precision and reasonable accuracy is 
possible for the brightest asteroseismic targets, with or with¬ 
out independent constraints on the radius and luminosity. It is 
also clear that our results are largely insensitive to fine details 
of the modeling strategy, but significant biases are possible 
when using simple power law corrections for surface effects 
and when neglecting diffusion and settling of heavy elements. 
Efforts are underway to define the character of these system¬ 
atic errors across a broad range of stellar properties, and to 
minimize or correct the model deficiencies that contribute to 
the biases. We look forward to additional opportunities to 
demonstrate the reliability of asteroseismic properties using 
observations of cluster members, binary systems, and inter¬ 
ferometric targets. 
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